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FOREWORD 


This  report  presents  methods  of  reducing  the  amount  of  computer 
storage  required  to  perform  the  computation  of  satellite  orbits. 

A  factor  of  six  reduction  will  be  realized  by  changing  current 
practice  at  the  Naval  Surface  Weapons  Center/Dahlgren  Laboratory  to 
that  suggested  here. 

The  study  was  performed  in  support  of  DoD  activities  utilizing 
near  earth  satellites. 


The  report  was  reviewed  by  R.  J.  Anderle,  Head,  Astronautics  and 
Geodesy  Division,  Naval  Surface  Weap.ns  Center/Dahlgren  Laboratory. 
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ABSTRACT 


This  study  was  undertaken  to  demonstrate  the  feasibility  of  storing 
just  the  position  segment  of  the  satellite  ephemeris  at  maximum  intervals 
of  time.  Position  of  the  satellite  at  times  not  on  the  ephemeris  would 
be  found  by  interpolation.  Satellite  velocity  and  acceleration  would 
be  calculated  via  numerical  differentiation  of  the  positions.  Partial 
derivatives  of  satellite  velocity  with  respect  to  orbit  parameters  would 
be  obtained  from  the  numerical  time  derivatives  of  the  position  partials 
with  respect  to  orbit  parameters.  Satellites  at  three  specific 
altitudes  were  considered:  (1)  a  low  altitude  satellite  (160  km), 

(2)  a  medium  altitude  (960  km),  (3)  a  high  altitude  satellite  (12000  km). 
To  achieve  acceptable  accuracy  standards,  the  optimum  interpolation 
order  and  storage  interval  as  a  function  of  satellite  altitude  are  as 
follows: 


Satellite 

Altitude 

Interpolation 

Interval 

Storage 

Interval 

low 

08 

90  sec 

med 

08 

180  sec 

high 

08 

900  sec 
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INTRODUCTION 

Given  the  initial  position  and  velocity  vector  of  a  satellite  at 
some  time  T,  a  Cowell  integration  routine  is  used  to  generate  position 
and  velocity  components  at  some  time  later  than  the  initial  time,  ihe 
integration  interval  used  is  a  function  of  the  satellite  altitude,  with 
low  altitude  satellites  being  integrated  at  smaller  intervals  because 
the  lower  satellites  are  more  sensitive  to  irregularities  m  the  earth  s 
gravitational  field.  See  Reference  1. 

Table  1  shows  how  the  Naval  Surface  Weapons  Center/Dahlgren  Laboratory 
stores  the  satellite  epnemeris.  Position,  velocity,  acceleration  and 
partial  derivatives  of  position  and  velocity  with  respect  to  ori 
parameters  are  given  at  each  integration  time  step.  Hie  satellite 
velocity  and  acceleration  are  time  derivatives  of  the  satellite 
position.  The  partial  derivatives  of  velocity  with  respect  to  orbi 
parameters  are  time  derivatives  of  the  respective  position  partials 
with  respect  to  orbit  parameters,  i.e. 


dx  _  d 


5x 

Tp- 


Thus  if  the  position  ephemeris  is  defined  to  be  that  part  of  the 
complete ’ephemeris  containing  satellite  positions  and  Partial  derivatives 
of  position  with  respect  to  orbit  parameters,  then  m  theory  it  is 
possible  to  reconstruct  the  complete  ephemeris  from  the  position 
eohemeris.  It  is  the  purpose  of  this  report  to  determine  exactly  this. 
Can  the  position  ephemeris,  stored  at  a-  Nerval  of  time  KT,  where  T 
is  the  integration  interval  and  K  is  the  positive  integer,  be  used  to 
reconstruct  the  complete  ephemeris  within  acceptable  (accuracy) 
tolerances? 
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The  partial  derivatives 
indicated  by  line  3 
(words  5  and  6)  and 
line  4  (words  1,2,3 
end  4)  are  repeated 
for  each  parameter 
to  be  altered 
(maximum  of  22 
parameters).  The 
number  of  parameters,  itr, 
appears  in  block  *1, 
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PROCEDURE 


Initial  conditions  for  three  representative  satellites  of  differing 
altitudes  were  integrated  forward  and  a  complete  ephemeris  was  generated 
for  each.  Two  trajectories  were  integrated  for  the  low  and  medium 
altitude  satellites.  The  low  altitude  satellite  at  15  and  30  seconds. 
The  medium  altitude  satellite  at  30  and  60  seconds.  The  integration 
interval  for  the  high  altitude  satellite  was  60  seconds.  lable  2 
shows  the  orbit  parameters  for  each  satel  Lite 

The  extended  storage  interval  KT,  is  obtained  by  sampling  every 
Kth  integration  time  line  from  the  nominal  ephemeris.  In  order  to 
reconstruct  the  ephemeris  to  its  original  denseness,  ie.  a  write  interval 
equal  to  the  integration  interval,  a  standard  centered  Lagrange 
interpolation  formula  of  the  form 

mfl 

f(t)  'K?-m  VS,fK  <l) 


is  used.  The  order,  M,  of the  interpolation  routine  is  equal  to  2m  +  2. 
Satellite  velocity  and  acceleration  are  obtained  from  the  position 
segment  of  the  ephemeris  by  differentiating  Equation  (1)  with  respect 
to  t  once  and  twice,  respectively.  See  Appendix  A  for  the  mathematical 
formulation.  In  similar  fashion  the  position  partials  are  reconstructed 
to  their  original  denseness  using  Equation  (1),  and  the  velocity  partials 
are  reconstructed  by  using  the  position  partials  and  the  time  derivative 
of  Equation  (1).  The  storage  interval  KT  and  the  interpolation  order 
M  were  varied  for  each  satellite  until  a  K  and  M  were  found  such  that 
the  complete  ephemeris  could  be  reproduced  from  the  position  segment 
to  within  the  following  tolerances  (in  any  one  component):  position; 

1.0X  lO^cm.  Velocity;  1.5  X  10"6  kn/'sec .  Acceleration;  1.5  X  10"  km/sec  . 

As  mentioned  above  the  initial  conditions  of  the  low  and  medium 
altitude  satellites  were  integrated  at  two  different  time  intervals. 

For  example,  consider  the  low  altitude  satellite  and  an  extended 
storage  interval  of  90  seconds.  If  the  trajectory  is  integrated  at  a 
30  second  time  step,  no  comparison  can  be  made  between  reconstructed 
ephemeris  values  and  the  actual  integrated  values  at  the  half  interval 
value  of  45  seconds.  Thus  the  necessity  of  having  a  15  second  integration 
interval.  On  the  other  hand  with  a  storage  interval  of  120  seconds, 
a  30  second  integration  interval  will  allow  comparison  at  half  step 
storage  intervals. 

Tables  3,  4,  5,  show  the  results  of  position,  velocity,  and 
acceleration  recovery  as  a  function  of  storage  interval  and  interpolation 
order,  '[he  Rq  column  shows  the  maximum  deviation  that  occurred  between 
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TABLE 
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SATELLITE 

PERIOD 

(min) 

INCLINATION 

(°) 

APOGEE  HEIGHT 
(nm) 

PERIGEE  HEIGHT 
(nm) 

Low 

88.66 

95.20 

147.02 

74.31 

Medium 

106.93 

89.25 

603.98 

563.08 

High 


469.14 


125.00 


7350.26 


7350.26 
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TABLE  4 

MEDIUM  ALTITUDE  SATELLITE 
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ORDER  __ 

INTERVAL 
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any  position  component  when  the  interpolated  position  was  differenced 
with  the  integrated  position.  Likewise  RD  and  RD  show  maximum  deviation 
for  velocity  and  acceleration  components  obtained  from  the  position  com¬ 
ponents.  Lastly  Rj  and  R-j-  show  the  maximum  deviation  in  any  one  com¬ 
ponent  between  interpolated  velocity  and  acceleration  obtained  using 
the  designated  storage  interval  and  the  value  of  velocity  and 
acceleration  from  the  nominal  complete  ephemeris.  All  maximum  deviations 
are  for  a  period  of  time  covering  more  than  one  revolution  of  the 
sa  tellite . 

Low  Altitude  Satellite 

The  optimum  interpolation  order  and  storage  interval  for  the  low 
altitude  satellite  is  (from  Table  3)  8^  order  and  90  seconds.  Accurate 
retrieval  cannot  be  made  using  a  120  second  storage  interval  even 
with  a  12t*1  order  routine.  Notice  also  that  in  almost  all  cases  the 
interpolated  velocities  and  accelerations  are  more  accurate  than  those 
obtained  by  numerical  differentiation  of  the  position  components. 

Medium  Altitude  Satellite 

The  optimum  case  for  the  medium  altitude  satellite  is  the  order 

routine  with  a  180  second  storage  interval.  Here  a  12th  order  routine 
and  a  240  second  interval  passes  the  required  tolerances.  However  it 
was  felt  that  the  computer  time  necessary  for  retrieval  using  this 
higher  order  outweighs  the  storage  advantage  of  the  longer  interval. 

High  Altitude  Satellite 

The  optimum  interpolation  order  and  storage  interval  for  the  high 
altitude  satellite  is  8^  order  and  900  seconds.  Unlike  the  case  of 
the  low  altitude  satellite,  the  high  altitude  satellite  shows  considerably 
more  accuracy  in  interpolated  velocity  and  acceleration  components  than 
their  corresponding  differentiated  counterparts. 

Appendix  B  contains  plots  of  the  resulting  differences  using 
various  storage  intervals  and  interpolated  orders  for  each  satellite. 

For  each  interpolation  order  and  storage  interval  there  is  a  set  of 
five  plots.  Interpolated  position  versus  nominal  position;  differentiated 
and  twice  differentiated  position  differences  giving  velocity  and 
acceleration  differences;  and  interpolated  velocity  and  acceleration 
versus  nominal  values.  Under  each  title  is  the  following  key 
(XX  0RD  YYY  SEC  ZZH) .  XX  gives  the  order  of  the  interpolation;  YYY 
the  storage  interval  and  ZZH  contains  the  integration  step  of  the  nominal 
trajectory  in  seconds  plus  the  altitude  of  the  satellite.  Thus  10  0RD 
180  SEC  30  M  would  identify  the  medium  altitude  satellite  with  the 
nominal  trajectory  integrated  at  30  seconds,  a  lO1"^  order  routine  and 
a  storage  interval  of  180  seconds. 
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The  next  area  of  consideration  involves  the  partial  derivatives. 

We  attempt  to  determine  to  what  degree  of  accuracy  the  partial  derivatives 
can  be  reconstructed  from  the  position  ephemeris  using  the  optimum 
interpolation  order  and  storage  interval  determined  above  for  each 
satellite.  Again  the  position  partials  are  reconstructed  at  the 
nominal  write  interval  using  Equation  (1),  and  the  velocity  partials 
are  reconstructed  using  the  position  partials  and  the  time  derivative 
of  Equation  (1). 

The  resultar  t  reconstructed  position  partial  derivatives  were 
normalized  and  compared  to  nominal  values  by  computing  the  following 
di f ference , 
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where  the  radical  subscripted  by  INTER  represents  interpolated  values 
and  the  entries  within  the  second  radical,  nominal  values.  The  velocity 
partials  were  compared  in  like  manner 
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Table  6  lists  the.  maximum  differences  between  reconstructed  normalized 
partials  and  their  normalized  nominal  counterparts.  3  UNITS,,  indicates 
a  maximum  deviation  of  3  units  in  the  eight  significant  place.  Maximum 
deviations  for  satellite  positions  ether  than  near  polar  or  equatorial 
regions  were  found  to  be  much  less  than  the  values  given  in  Table  6. 

In  an  attempt  to  give  physical  significance  to  the  entries  of  Table  6, 
the  foilwing  test  was  conducted.  In  the  case  of  the  low  altitude  satel¬ 
lite  two  propogated  trajectories  were  constructed.  (A  propogated 
trajectory  is  constructed  using  the  partial  derivatives  of  position  and 
velocity  with  respect  to  orbit  parameters  and  a  set  of  improvements  to 
the  orbital  elements  to  make  corrections  to  the  original  trajectory.) 

The  first  propogated  trajectory  was  formed  using  the  nominal  position 
ephemeris  written  at  a  15  second  interval  and  a  representative  set  of 
changes  in  orbital  parameters.  The  second  propogated  trajectory  for 
the  low  altitude  satellite  was  formed  using  the  nominal  low  altitude 
ephemeris  written  at  90  seconds  plus  the  same  se-t  of  representative 
changes  to  orbital  parameters.  An  order  Lagrangian  interpolation 
routine  was  used  on  this  trajectory  to  write  position  partials  every 
15  seconds.  Figure  1  shows  the  total  position  deviation  between  these 
two  propogated  trajectories.  The  maximum  difference  is  4  cm. 


In  like  mai  ner  two  propogated  trajectories  were  created  for  the 
medium  altitude  satellite.  The  first  was  computed  using  the  nominal 
medium  altitude  trajectory  written  at  60  seconds  and  a  representative 
set  of  changes  to  orbital  parameters.  The  second  propogated  trajectory 
used  the  nominal  medium  altitude  trajectory  written  at  180  seconds, 
figure  2  shows  the  total  difference  between  these  two  trajectories. 
Maximum  deviation  is  below  the  1  cm  level. 

In  each  case  the  difference  between  the  two  trajectories  is  small 
adding  redence  to  the  fact  that  the  interpolation  orders  and  storage 
intervals  listed  above  are  indeed  optimum  not  only  for  position 
velocity  and  acceleration  but  also  for  partial  derivatives. 

As  a  result  of  the  above  findings  it  is  now  possible  to  store  iust 
the  position  segment  of  the  ephemeris  (i.e.  satellite  position  and 
partial  derivatives  of  position  wi’th  respect  to  dynamic  parameters 
at  a  storage  interval  greater  than  the  integration  interval.  This 
results  m  a  tremendous  savings  in  storage  requirements.  Omitting 
velocity  components  and  partial  derivatives  of  velocity  with  respect  to 
ynamic  parameters  results  in  storage  requirements  being  reduced  by  a 
factor  of  2.  In  the  case  of  the  low  and  medium  altitude  satellites' 
ephemeris  presently  stored  at  30  seconds  and  60  seconds,  respectively 
scoring  the  position  ephemeris  at  the  optimum  storage  intervals  results 
in  total  storage  requirements  being  reduced  by  a  factor  of  6.  Even 
greater  savings  result  in  the  case  of  the  high  altitude  satellite. 
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APPENDIX  A 


MATHEMATICAL  FORMULATION 


Ma thematical  Formulation 


llIMfilumin. 


The  standard  centered  Lagrange  interpolation  formula  is: 

m+1 

(1) 


where: 
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At 


To  avoid  loss  of  accuracy  due  to  a  zero  divisor  as  s-0  or  s~*l , 
Equation  1  is  rewritten  as: 
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where  (s+m)n^,  (s+m) , . . (s+1) (s-2) (s-3) . . . (s-m-1)  . 
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For  the  1st  derivative  we  differentiate  (2)  and  obtain 


m+1  p 

*  „  1  v  m,K 

£(t)  At  '  K?-m  (s-K) 

I  K^O ,  1 


S  - 


(s+m) 


[2m+2] 


(s-K) 


J 


[2m+2] 


m+1 


+  <^)0;rJ  Pm,0£o(lf(S-1)P)+  Pm,lfl(1+sn 

where  S-  (s+m)  2”t2]P  +  (s+mj“2]  (2s-l)  and  P  -  2.^ 

K^0,1 


(3) 


A-l 


. . . . . 


s-m- 1)  , 


(2) 


The  second  derivative  is  found  by  differentiating  (3).  The  result  is 
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APPENDIX  B 

PLOTS  OF  RESULTING  DIFFERENCES 
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